Cosmological particle creation in the lab? 
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2 2 SCATTERING ANALOGY 

1 Introduction 

One of the most striking examples for the production of particles out of the quantum 
vacuum due to external conditions is cosmological particle creation, which is caused by 
the expansion or contraction of the Universe. Already in 1939, Schrodinger understood 
that the cosmic evolution could lead to a mixing of positive and negative frequencies and 
that this "would mean production or annihilation of matter, merely by the expansion" 
[Schrodinger, 1939]. Later this phenomenon was derived via more modern techniques of 
quantum field theory in curved space-times by Parker [Parker, 1968] (who apparently 
was not aware of Schrodinger's work) and subsequently has been studied in numerous 
publications, see, e.g., [Birrell & Davies, 1982; Fulling, 1989; Wald, 1994]. Even though 
cosmological particle creation typically occurs on extremely large length scales, it is one 
of the very few examples for such fundamental effects where we actually may have ob- 
servational evidence: According to the inflationary model of cosmology, the seeds for 
the anisotropies in the cosmic microwave background (CMB) and basically all large scale 
structures stem from this effect, see Section In this Chapter, we shall provide a brief 
discussion of this phenomenon and sketch a possibility for an experimental realization via 
an analogue in the laboratory. 



2 Scattering analogy 

For simplicity, let us consider a massive scalar field $ in the 1+1 dimensional Friedmann- 
Robert son- Walker metric with scale factor a{r) 

ds 2 = dr 2 - a 2 {r) dx 2 = a 2 ( V ) [drf - dx 2 ] , (1) 

where r is the proper (co-moving) time and r] the conformal time. The latter co-ordinate 
is more convenient for our purpose since the wave equation simplifies to 

d 2 d 2 



drj 2 dx 2 



a {rj)m $(t),x) . (2) 



In the massless case m — 0, the scalar field is conformally invariant (in 1+1 dimensions) 
and thus the expansion does only create particles for m > 0. After a spatial Fourier 
transform, we find that each mode 4>k(v) behaves like a harmonic oscillator with a time- 
dependent potential 

|^ + fi 2 (t)W) = 0, (3) 



with k 2 + a 2 (r])m 2 — > Q 2 (t) and r] —$■ t. There is yet another analogy which might 
be interesting to notice. If we compare the above equation to a Schrodinger scattering 
problem in one spatial dimension 

■T- T~2 + V ( x )) *(*) = m ^ . ( 4 ) 

we find that is has precisely the same form after identifying t x, 0(t) <h> ^(x), and 
f2 2 (£) -H- 2m[E — V(x)]. Note that f2 2 is always greater than zero in our case - which 
corresponds to propagation over the barrier E > V{x). If Vt 2 were less than zero over some 
region in time, one would have a barrier penetration (i.e., tunnelling) problem E < V(x). 
With the condition that in the past the field has the form e lQint , in the future the solution 
would be ae' tn ° utt + /3e~ in ° utt due to scattering from the region where Q 2 < 0. This would 
correspond to particle creation with probability proportional to |/3| 2 . However even if 
Q 2 > everywhere there will still be some scattering (above the barrier). 
In order to derive the cosmological particle creation, we can study a positive pseudo-norm 
solution of Eq. ([3]) which initially behaves as e~ in ' int and finally evolves into a mixture of 
positive and negative pseudo-norm solutions - which is in this case equivalent to positive 
and negative frequencies ae~ lUoutt +/3e + outt (assuming that Q is constant asymptotically). 
In the Schrodinger scattering problem, the initial solution e~ lQt could be identified with a 
left-moving wave on the left-hand side of the potential "barrier" while the final solution 
ae~ tnt +(3e +int would then correspond to a mixture of left-moving ae~ lQt and right-moving 
(3e +tnt waves on the right- hand- side. As a consequence, the Bogoliubov coefficients a and 
f3 are related to the reflection R and transmission T coefficients via a = 1/T and /3 = R/T. 
In this way, the Bogoliubov relation \a\ 2 — |/3| 2 = 1 is equivalent to the conservation law 
\R\ 2 + \T\ 2 = 1 for the Schrodinger scattering problem. The probability for particle 
creation can be inferred from the expectation value of the number of final particles in the 
initial vacuum state which reads (0 in | h out |0 in ) = \(3\ 2 . 

3 WKB analysis 

In order to actually calculate or estimate the Bogoliubov coefficients, let us re- write Eq. (J3]) 
in a first-order form via introducing the phase-space vector u and the matrix M 

= *=( ° 2 1 )-(*)=M-u. (5) 

^ -ff (t) o ) \ + ) yl 

If we define an inner product via 

(u\u') = i(^2 n i — n l n 2) > (6) 




4 3 WKB ANALYSIS 

we find that the inner product of two solutions u and u' of Eq. (jSj) is conserved 

The split of a solution into positive and negative frequencies (i.e., positive and negative 
pseudo-norm) corresponds to a decomposition in the instantaneous eigen-basis of the 
matrix 

M u± = ±iOu± . (8) 

Choosing the usual normalization u± = (1, ±iO) T /V20, we find 

(u+|ii+) = l, (u-\u-) — — 1, (u+|it_) = . (9) 

At each time t, we may expand a given solution u(t) of Eq. (jSJ) into the instantaneous 
eigen- vectors 

ti(t) = a(t)e iv ®u+{t) + /3(t)e-^ (t) u_(t) , (10) 

where the pre-factors are now defined as time-dependent Bogoliubov coefficients a(t) and 
0(t). It is useful to separate out the oscillatory part with the WKB phase 
t 
(p(t) = f dt'tt(t'). (11) 

— oo 

Now we may insert the expansion f lTOj) into the equation of motion (JSJ) and project it with 
the inner product (jSJ) onto the eigen-vectors u± which gives 

* = ^-*^ = £«"'«■ < 12 > 

due to (ti_|tit + ) = (l/(2Q,) and (tt + |tt_) = —0/(20) while (tt + |tt + ) = (tt_|tt_) = 0. 
This equation ( 1T2|) is still exact and very hard to solve analytically - except in very special 
cases. It can be solved formally by a iterative integral equation 

t 
«„ + i = a m + J dt'^e-^p n (t'), 

— oo 

t 
/?n + i = (3 m + Jdt'^e-^a n (t'). (13) 

— oo 

It can be shown that this iteration converges to the exact solution for well-behaved 0(t) 
[Braid, 1970]. Standard perturbation theory would then correspond to cutting off this 



iteration at a finite order, which can be justified if Q(t) changes only very little. For 
the scalar field in Eq. fl2]) this perturbative treatment should be applicable in the ultra- 
relativistic limit, i.e., as long as the mass is much smaller than the wave-number. 
In many cases, however, another approximation - the WKB method - is more useful. This 
method can be applied if the rate of change of £l(i), e.g., the expansion of the universe, 
is much slower than the internal frequency Q(t) itself. Writing 

fi(t) = fio/M) , (14) 

with some dimensionless function / of order one, the WKB limit corresponds to Qq ^> u. 
In terms of the reflection coefficient R = /3/a mentioned earlier, we get 

R=-^r (e 2i * - R 2 e~ 2 ^) , (15) 

which is known as Riccati equation. Again, this equation is still exact but unfortunately 
non-linear. Neglecting the quadratic term R 2 would bring us back to perturbation theory. 
In the WKB-limit, the phase factors e ±2l(p are rapidly oscillating and the magnitude of R 
can be estimated by going to the complex plane. Re-writing the Riccati equation (115]) as 

^ = 1 (e 2 - - R 2 e~ 2 ^ ^ , (16) 

dip 2 v ; dp ' K ' 

we may use an analytic continuation p — > <p + i\ to see that R becomes exponentially sup- 
pressed R ~ e~ 2x . How strongly it is suppressed depends on the point where the analytic 
continuation breaks down. Since e ±2tip is analytic everywhere, this will be determined by 
the term lnfi. Typically, the first non-analytic points £* encountered are the zeros of Q, 
i.e., where Cl(t*) = 0. In the case of barrier reflection, these points where Q = 0, i.e., 
where V = E, lie on the real axis and correspond to the classical turning points in WKB. 
In our case, we have scattering above the barrier and thus these points become complex 
- but are still analogous to the classical turning points in WKB. Consequently, we find_J 



a 



R = - ~ e ~ 2x * = exp { -23 ( / dt'tt(t') ) \ . (17) 



If there is more than one turning point, the one with the smallest x* > 0, i.e., closest to 
the real axis (in the complex <£>-plane) dominates. If these multiple turning points have 
similar x* > 0, there can be interference effects between the different contributions, see, 
e.g., [Dumlu & Dunne, 2010]. 



^n fact, it can be shown that Eq. (fl7|) becomes exact in the adiabatic limit u/Cl J, 0, i.e., the pre-factor 
in front of the exponent tends to one, see, e.g., [Massar & Parentani, 1998; Davis & Pechukas, 1976]. 



6 4 ADIABATIC EXPANSION AND ITS BREAKDOWN 

4 Adiabatic expansion and its breakdown 

Note that we could repeat steps (jSJ) till (|12p and expand the solution u(t) into the first- 
order adiabatic eigen-states instead of the instantaneous eigen- vectors u±. To this end, 
let us re-write (fT2l) as 



P(t)e-*® J ' ' I fi/(2fi) -ifi J ' I /3(t)e~^W 

The eigen-vectors of the matrix N are the first-order adiabatic eigen-states w± and the 
eigen-frequencies N ■ w± = ±if2 a dio± are renormalized to 




Cl 2 
1_ 4^- (19) 

Assuming a- m = 1 and f3- m = 0, the system stays in the adiabatic eigen-state w + to lowest 
order in u/Q and we get 

<*) = 1 + O (g) , P{t) = -1| + O g) . ,20) 

This adiabatic expansion into powers of uj/Vt can be continued and gives terms like Cl 2 /Q A 
and Cl/Q 3 to the next order in uj/Vt (see below). One should stress that this expansion is 
not the same as in ( Tl3|) since it is local - i.e., only contains time-derivatives - while ( TT31) 
is global - i.e., contains time- integrals. Since all terms of the adiabatic expansion ( 1201) 
are local, they cannot describe particle creation - which depends on the whole history 
of Q(t). In terms of the adiabatic expansion into powers of ou/Q, particle creation is a 
non-perturbative effect, i.e., it is exponentially suppressed, see Eq. ( TT7|) 

E-expj-O^U , (21) 

and thus cannot be found be a Taylor expansion into powers of lj/Q. For any finite ratio 
of u/Q, this also means that the adiabatic expansion (into powers of oj/Q) must break 
down at some point. To make this argument more precise, let us re- write Eq. f fl8|) in yet 
another form 

^ = N . w = A ( i ™ h W s ' nh < 2 « ).„. (22) 

dt \ sinh(20 -zcosh(2£) J 

In this representation, the eigen- values of N are given by ±z'A and the eigen-vectors read 

( COSn £ \ ( Sinh ^ \ /OQ\ 

w += \ . . , , > w- = . , , ■ ( 23 ) 





Decomposing the solution w(t) into these eigen- vectors 

w(t) = a(t)w+(t) + b(t)w_(t) , (24) 

and using w + = £«;„ as well as w = £w + , we find 

(25) 

This is the same form as Eq. ( 12 2 j) if we change A and £ accordingly. Thus, by repeating 
this procedure, we get the iteration law 







A„+i = V Al - H , Cn+i = -r arctanh — . (26) 



By this iteration, we go higher and higher up in the adiabatic expansion since £ n always 
acquires an additional factor of u/Q. Thus, for u <C Q, the values of C, n quickly decay 
with a power-law £„ = 0([u;/f2] n ) initially. As we go up in this expansion, however, the 
effective rate of change of £ n increases. For example, if Q(t) has one global maximum (or 
minimum) and otherwise no structure, the time-derivative Cl/(2Q 2 ) = tanh(2£i) has two 
extremal points and a zero in between. By taking higher and higher time derivatives, more 
and more extremal points and a zeros arise and thus the effective frequency u^ of £ n (£) 
increases roughly linearly with the number n of iterations uf = 0{nuS). Furthermore, 
the adiabatically renormalized eigen-values A n decrease with each iteration. Thus, after 
approximately n = 0(Q/oj) iterations, the effective frequency u^ becomes comparable to 
the internal frequency A n . At that point, the adiabatic expansion starts to break down. 
Estimating the order of magnitude of £ n at that order gives 



6> = °([n]) = °([n] ) 



(27) 



Since the effective external wf and internal A n frequencies are comparable and £ n is 
very small, we may just use perturbation theory to estimate j3 and we get j3 = 0(£ n ), 
i.e., the same exponential suppression as in Eq. (I2"T1) . If we would continue the iteration 
beyond that order, the £ n would start to increase again - which the usual situation in an 
asymptotic expansion, see Figure [U Carrying on the iteration too far beyond this point, 
the ^ exceed the A^ and thus we have barrier penetration instead of propagation over 
the barrier (as occurs for all orders below this value of n). In this procedure, it is this 
barrier penetration which gives the mixing of positive and negative pseudo-norm, and the 
creation of particles. Were the system to remain as propagation over the barrier for all 
orders n in this adiabatic expansion, one would have no particle creation. 
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Iteration Number 



Figure 1: Sketch of the effective external frequencies cu° ff (crosses) and amplitudes £ n (solid 
line) depending on the iteration number n obtained numerically for a concrete example. 
One can observe that uf grows approximately linearly with n while £ n first decreases but 
later (for n > 5) increases again. 



5 Example: inflation 

As an illustrative example, let us consider a minimally coupled massive scalar field in 
3+1 dimensions - which could be the inflaton field (according to our standard model of 
cosmology). Again, we start with the Friedmann- Robert son- Walker metric ([T]) with a 
scale factor a{r) and obtain the equation of motion 



J_ d_ 3 d_ 
aHr) dr ° [T) 8r 



a 2 (r) 



m 



$ = 0. 



m 
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Rescaling the field 0(r, r) = l5(r)&(T,r) with U(t) = a 3 ^ 2 (r) and applying a spatial 
Fourier transform, we obtain the same form as in Eq. ([3]) 

d 2 k 2 2 1 d 2 U(r)\ , n 

In the standard scenario of inflation, the space-time can be described by the de Sitter 
metric a(r) = exp{H r} to a very good approximation, where H is the Hubble parameter. 
In this case, the effective potential 15/15 just becomes a constant (3H/2) 2 and the frequency 
f2(r) reads 

« 2 M^W-^. (30) 

Inserting a{r) = exp{Hr}, we see that modes with different k- values follow the same 
evolution - just translated in time. (This fact is related to the scale invariance of the 
created k spectrum.) Initially, this frequency is dominated by the k term and we have 
Ct/ft = —H which means that we are in the WKB regime Q/Q <C Q. However, due to 
the cosmological red-shift, this k term decreases with time until the other terms become 
relevant. Then the behavior of the modes depends on the ratio m/H. For m ~^> H, 
the modes remain adiabatic (i.e., stay in the WKB regime) and thus particle creation is 
exponentially suppressed. If m and H are not very different, but still m > 3H/2 holds, 
the modes are adiabatic again for large times - but for intermediated times, the WKB 
expansion breaks down, leading to a moderate particle creation. For m < 3H/2, on the 
other hand - which is (or was) supposed to be the case during inflation - the frequency 
fi(r) goes to zero at some time and becomes imaginary afterwards. This means that we 
get a barrier penetration (tunneling) problem where the modes (pk{ T ) do not oscillate but 
evolve exponentially in time 4>k{ T ) oc exp{±r ^/9H 2 / A — m 2 }. Here one should remember 
that the original field does not grow exponentially due to the re-scaling with the additional 
factor 15(t) = a?' 2 {j). This behavior persists until the barrier vanishes, i.e., the expansion 
slows down (at the end of the inflationary period) and thus the effective potential 15/15 
drops below the mass term. After that, the modes start oscillating again. However, in 
view of the barrier penetration (tunneling) over a relatively long time (distance), we get 
reflection coefficients R which are not small but extremely close to unity R « 1. This 
means that the Bogoliubov coefficients a and (3 are huge - i.e., that we have created 
a tremendous amount of particles out of the initial vacuum fluctuations. According to 
our understanding, precisely this effect is responsible for the creation of the seeds for all 
structures in our Universe. Perhaps the most direct signatures of this effect are still visible 
today in the anisotropics of the cosmic microwave background radiation. 
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An alternative picture of the mode evolution in terms of a damped harmonic oscillator 
can be obtained from the original field in Eq. ( 128]) 

Initially, the term e~ 2HT k 2 dominates and the modes oscillate. Assuming m <^ H (which 
is related to the slow-roll condition of inflation), the damping term dominates for late 
times and we get a strongly over-damped oscillator, whose dynamics is basically frozen 
(like a pendulum in a very sticky liquid). The transition happens when H ~ ke~ Hr , i.e., 
when the physical wavelength A = 2ire HT /k exceeds the de Sitter horizon oc 1/H due to 
the cosmo logical expansion e Hr . After that, crest and trough of a wave lose causal contact 
and cannot exchange energy any more - that's why the oscillations effectively stops. 
As a final remark, we stress that this enormous particle creation effect is facilitated by 
the rapid (here: exponential) expansion and the resulting stretching of wavelengths over 
many many orders of magnitude (i.e., the extremely large red-shift). Therefore, a final 
mode with a moderate wavelength originated from waves with extremely short wave- 
lengths initially. Formally, these initial wavelengths could be easily far shorter than the 
Planck length. However, on these scales one would expect deviations from the theory 
of quantum fields in classical space-times we used to derive these effects. On the other 
hand, this problem is not only negative - it might open up the possibility to actually 
see signatures of new (Planckian) physics in high-precision measurements of the cosmic 
microwave background radiation, for example. 



6 Laboratory analogues 

Apart from the observation evidence in the anisotropics of the cosmic microwave back- 
ground radiation mentioned above, one may study the phenomenon of cosmological par- 
ticle creation experimentally by means of suitable laboratory analogues, see, e.g., [Unruh, 
1981; Barcelo, Liberati, & Visser, 2011]. The are two major possibilities to mimic the 
expansion or contraction of the Universe - a medium at rest with time-dependent proper- 
ties (such as the propagation speed of the quasi-particles) or an expanding medium. Let 
us start with the former option and consider linearized and scalar quasi-particles (e.g., 
sound waves) with low energies and momenta propagating in a spatially homogeneous and 
isotropic medium. Under these conditions, their dynamics is governed by the low-energy 
effective action 

£ cff = \ (a 2 (t)<p 2 + b 2 {t)<p 2 + c 2 (t)[V0] 2 ) + 0{^) + 0(d 3 ) . (32) 
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Here we assume positive a 2 and non-negative b 2 and c 2 for stability. The factor a 2 (t) can 
be eliminated by suitable re-scaling of the time co-ordinate. Then, after a spatial Fourier 
transform, we obtain the same form as in Eq. (|2|). The quasi-particle excitations <fi in 
such a medium behave in the same way as a scalar field in an expanding or contracting 
Universe with a possibly time-dependent potential (mass) term oc b 2 (t)(f) 2 . In order to 
avoid this additional time- dependence of the potential (mass) term, the factors b and c 
must obey special conditions. For example, Goldstone modes with 6 = correspond to a 
massless scalar field in 3+1 dimensions - whereas the case of constant c is analogous to a 
massive scalar field in 1+1 dimensions. 

As one would intuitively expect, the expansion or contraction of the Universe can also be 
mimicked by an expanding or contracting medium. Due to local Galilee invariance, such 
a medium can also be effectively spatially homogeneous and isotropic as in Eq. (13"2"j) when 
described in terms of co-moving co-ordinates. For a quite detailed list of references, see 
[Barcelo, Liberati, & Visser, 2011]. 

There are basically three major experimental challenges for observing the analogue of 
cosmological particle creation in the laboratory. First, the initial temperature should be 
low enough such that the particles are produced due to quantum rather than thermal 
fluctuations. Second, one must be able to generate a time- dependence (e.g., expansion of 
the medium) during which the effective action in Eq. (132]) remains valid (in some sense) 
but which is also sufficiently rapid to create particles. Third, one must be able to detect 
the created particles and to distinguish them from the radiation stemming from other 
sources. For trapped ions, for example (see, e.g., [Schiitzhold et al, 2007]), the first and 
third point (i.e., cooling and detection) is experimental state of the art, while a sufficiently 
rapid but still controlled expansion/contraction of the ion trap presents difficulties. For 
Bose-Einstein condensates (see, e.g., [Barcelo, Liberati, & Visser, 2011] and references 
therein), on the other hand, the first and third points are the main obstacles. 
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